Quantum dynamical phase transition in a system with many-body interactions 
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We introduce a microscopic Hamiltonian model of a two level system with many-body interactions 
with an environment whose excitation dynamics is fully solved within the Keldysh formalism. If 
a particle starts in one of the states of the isolated system, the return probability oscillates with 
the Rabi frequency ujq- For weak interactions with the environment 1/t se < 2uiq, we find a slower 
oscillation whose amplitude decays with a rate l/rj, = l/(2r SE ). However, beyond a finite critical 
interaction with the environment, l/r SE > 2wo, the decay rate becomes l/r^ oc ljq t se . The 
oscillation period diverges showing a quantum dynamical phase transition to a Quantum Zeno phase. 



Ideal quantum information processing (QIP) involves 
manipulating a system's Hamiltonian. In practice, the 
interactions with an environment^ perturb the evolution, 
smoothly degrading the quantum interferences within a 
"decoherence" rate, 1/ta. Although one expects l/r^ to 
be proportional to the system-environment (SE) interac- 
tion rate 1/tse, there are conditions where l/r^ does not 
dependent on h— . This phenomenon was interpreted^^ 
as the onset of a Lyapunov phase, where the decay rate is 
the Lyapunov exponent A characterizing the complexity 
of the classical system. The description of this transition, 
1/t0 = min [1/tse, A], requires evaluation of the observ- 
ables beyond perturbation theory^. We will show that 
a dynamical transition also occurs in a swapping gate: 
a system that jumps between two equivalent states, A 
and B, when the coupling Vab is turned on£i Start- 
ing on state A, the return probability oscillates with the 
Rabi frequency u>q — 2V^b/?L Since each state inter- 
acts with the environment at a rate 1/t se , weak inter- 
actions (1/t se < 2uo) produce a slightly slower oscilla- 
tion which decays at a rate I/t^ = 1/(2t se ). However, 
the swapping frequency is non- analytic on the interac- 
tion rate and at a critical strength l/rg E = 2ojq, the 
oscillation freezes indicating a transition to a new dy- 
namical regime. The initial state now decays to equi- 
librium at a slow rate 1/ T x w o t se which cancels for 
strong SE interaction. This last regime can be seen as 
a Quantum Zeno phase, where the dynamics is inhibited 
by frequent "observations'— by the environment. Such 
quantum freeze can arise as a pure dynamical process 
governed by strictly unitary evolutions2iS. 

We consider a "system" with a single electron occu- 
pying one of two coupled states, A or B, each interact- 
ing with a corresponding electron reservoir (the "envi- 
ronment"). The total system, represented in Fig. Q] a) 
has the Hamiltonian H = Hs+'He+'Hse, where the first 
term is 



= E A c\c A + E B c%c B - V AB {c\c b + c+c A ) . (1) 



Here, c±(cA are the standard fcrmionic creation (destruc- 
tion) operators. -E 1 , stands for the energy of the i-th local 
state whose spin index is omitted. Vab yields the natural 



frequency, luq = 2Vab / h. The environment is 
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where the sums on negative (positive) index describe a 
semi-infinite chain to the left (right) acting as a reservoir. 
E—i = El and E\ = Er are site energies while V-i=Vl 
and V\ = Vr are adjacent site hoppings. The system- 
environment interaction is modeled with a through-space 
interaction 
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The first line represents the Coulomb interaction of an 
electron in state B with an electron in the first site of 
reservoir to the right, [/g^is the standard direct inte- 
gral and J7grt ^ s the exchange one. Analogously, the 
second line is the interaction with the reservoir to the 
left. 

A complete norm preserving solution requires the eval- 
uation of the reduced particle and hole density functions 
G§(fe,ti) = i(*|ct(t x ) c,^)]*) and G% (f 2 ,ti) = 
— H (*l £j (*i) I*) that describe temporal and spa- 

cial correlations. Here, the creation and destruction op- 
erators are in the Heisenberg representation and = 
c\c B | "Jo) is an initial non-equilibrium many-body state 
built from the non-interacting equilibrium |^o)- The re- 
tarded Green's function G^(t 2 ,t\) = [Gjj (ti, h)] = 
(i a , k) [G> (tz, h)- G< (t 2 , tx)] describes the probabil- 
ity amplitude of finding an electron at site i after placing 
it at site j and letting it evolve under the total Hamil- 
tonian for a time t% — t\. By restricting the analysis to 
i,j G {A,B}, G < (t,t) is the single particle 2x2 den- 
sity matrix and G R (£2^1) is an effective evolution op- 
erator in this reduced space. In absence of SE interac- 
tion, the Green's function is easily evaluated in its en- 
ergy representation G 0R (e) = J G 0R (t) exp[iet / h]dt = 
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FIG. 1: (Color online) a) System-environment scheme. 
Squares and circles represent system and environment states 
respectively. Solid lines are tunneling matrix elements while 
dashed lines are through-space Coulomb interactions, b) The 
relevant self-enegy diagrams in a local basis. Lines are ex- 
act Green's functions of the system and the environment in 
absence of Coulomb interactions (dashed lines), c) Retarded 
Green's function of the system-environment at site B. It con- 
tains the Coulomb interaction through the self-energy correc- 
tion of b) to infinite order, d) Particle density function at site 
B. The dashed lines represent local interactions in time and 
space. 



[el — Hs] 1 . Conversely, the interacting Green's function 
defines the reduced effective Hamiltonian and the self- 
energies S R (e)ii, Heff.(e) = el [G R (e)] _1 = H s + 
S R (e), where the exact perturbed dynamics is contained 
in the nonlinear dependence of the self-energies £ R on e. 
For infinite reservoirs Re E R (e°) represents the "shift" 
of the system's eigen-energies e° v and — Im£ R (e°) /K = 
l/(2r gE ) accounts for their "decay rate" into collective 
SE eigenstates in agreement with a self-consistent Fermi 
Golden Rule (FGR)i£, i.e. the evolution with H e g. is 
non-unitary. 

The complete dynamics will be obtained resorting to 
the Keldysh formalismi&ii. This allows the evaluation 
of the relevant density-density correlations within a norm 
conserving scheme^: 

G< (t 2 , h) = h 2 G R (t 2 , 0) G< (0, 0) G A (0, t x ) + 

r r di fc d^c R (t 2) t k ) s< (t k , u) g a (t h h) . (4) 

Jo Jo 

The first term stands for the "coherent" evolution 
while second term contains "incoherent reinjections" , de- 
scribed by the injection self-energy, S < , that compen- 
sates any leak from the coherent evolution—. Solv- 
ing Eq. J3J requires the particle (hole) injection and 
retarded self-energies, (ii, i 2 ) and Y> R (t\,t 2 ) = 



6{ti,t 2 ) [£> (t 2 ,h) - £< (t 2 ,ti)]. For this, we use a per- 
turbative expansion on TCse- The first order gives the 
standard Hartree-Fock energy corrections which, being 
real, do not contribute to S < . The second order termsi^, 
sketched in Fig. \J\b), contribute to E^, and in the space 
representation: 



h 2 



= \U i fGf s (t k ,t l )Gf s (t l ,t k )G? i (t k ,t l )6 ij , 

(5) 

where (i,s) G {(A, L) , (B, R)} and hence s stands for 
the surface site of the environment. The net interaction 
between an electron in the system and one in the reser- 
voir is Uis =— 2[/£ dlr '' ) + [/£° xch '\where the direct term 
contributes with a fermion loop with an extra spin sum- 
mation. Notice that self-energy diagrams shown in Fig. 
n b) describe an electron exciting an electron- hole pair 
in the environment and later destroying it. The evalu- 
ation of these processes requires accounting for the dif- 
ferent time scales for the propagation of excitations in 
the system and reservoirs. We resort to time-energy 
variables^: U = k t , the physical time, and Sti = 
tk ~ tli which characterizes the quantum correlations. 
The integrand in Eq. Q), when t 2 = t\ = t becomes 
G R (t, U + SU/2) £< (8U, U) G A (U - Sti/2, t). Since SU 
is related to the energy e through a Fourier transform 
(FT), in equilibrium, 



G< s (e,U)=i2n N s (s) k(e,*i)> 
G> (e,ti) = -i2n N s (e) [l-&(e,t< 



(6) 



Here N s (e) is the local density of states (LDoS) at the 
surface of the reservoir and f s (e, ti) = | is the occupa- 
tion factor in the high temperature limit (ksT 3> V s ). 



Replacing the LDoSi^ N s (e) 
Eq 



l/(nV s . 



1 - 



2\:. 



©, and doing the FT— one gets Gf tS (Sti,ti) = 

' V ' Jl St ^ ^ i ' wnere Ji i s the Bessel function of first 
order. Replacing in Eq. (JHJ 
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| Ji [2V s 6ti/H] I (V s Sti)\ decays in a time scale H/V s 
which, in the wide-band (WB) limit V s » Vab, is 
much shorter than H/Vab, the time scale of GfASti,ti). 
Hence, the main contribution to the integral on Sti in 
Eq. (JUl is obtained replacing Gf l (St i ,t l ) by G^(0,U). 
The same consideration holds for G R (t, U + Sti/2) and 
G A (ti — Sti/2, t) which are replaced by G R (i,ii) and 
G A (ti,t). Then, the dependence on Sti enters only 

through Eg (SU,U) yielding 



(ti 



\Ui, 



6tt V s 



(8) 
(9) 
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Hence, in the WB limit, the Keldysh self-energy of Eq. 

becomes local in space and time and no further struc- 
ture from higher order terms is admisible. This is repre- 
sented as a collapse of sucessive pairs of black dots in Fig. 
Q b) into a single point. The expansion represented in 
Fig. \i\c) an d that of Fig. \i\d) become exact in this limit. 

We assume Ea = Eb = -El = Er = and the sym- 
metry condition \U AL \ 2 /V L = \U B r\ 2 /Vr. From Eq. @ 
we obtain the decay rates 



where a 2 , 
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coinciding with the FGR. The WB limit can be relaxed 
because the FGR holds when time t is in the range 12 
ts < t < tn ~ ah/TsE ^(V^/rss)) where a depends 
on the van Hove singularities of the spectral density 
J s (e) = f N s (e - e')N s (e')de' . Here, t s = hj s (0) ~ 
fo/V^ is the survival time of an electron-hole excita- 
tion at the surface site and tji characterizes the transi- 
tion to a power law decay arrising from memory effects. 
Hence, as long as Tse, Vab <C V s , the FGR is valid for 
times much longer than H/Tse- Under these conditions, 
H e ff. does not depend on e and the propagator has a 
simple dependence on t as G R (t) = G 0R (t) e ~ r ^ t / r \ 
where, G° A R A (t) = G™(t) = \ cos (ft) and Gf B (t) = 
G Q B R A (t) = isin(^ft). Eq. Qi becomes, 



G< (i, t) = h 2 G 0R (t) G< (0, 0) G 0A (-t) e-*/ TSE + 

/ di;G 0R (t - U) S< (U) G 0A (ti - i) e -(*-*i)A*» 
Jo 



(11) 



a generalized Landauer-Biittiker equatio: 



,15,17 



The 



tial condition has the state A occupied: jGf^ (0,0) = 
SiA^Aj- Replacing Eq. |j5Jl into Eq. iJHJ, identifying the 
interaction rate of Eq. (fTU|) we get two coupled equations 
for G AA and G BB 



hG 
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(12) 



The first term is the probability that a particle initially at 
site ^4 is found in site A (or B) at time i having survived 
the interactions with the environment. The second and 
third terms describe particles whose last interaction with 
the environment, at time ti, occured at site A and B 
respectively. The Laplace transform yields: 



}(*AA (*i *) = l+ a o cos [{u + in)t 



-*/(2 r SE ) 



(13) 



Noticeably, in the first term of Eq. (|ll|l the environment, 
though giving the exponential decay, does not affect the 
frequency. Modification of lu requires the dynamical feed- 
back. 

The effect of lateral chains on the two state system can 
produce observables with non-linear dependences on -ffsE 
which could account for a cross over among the limiting 
dynamical regimes. However, we find a non-analyticity in 
these functions enabled by the infinite degrees of freedom 
of the environment^ (i.e. the thermodynamic limit). 
Here, they are incorporated through the imaginary part 
of the self-energy, H/tse, i.e. the FGR. Hence, the non- 
analyticity of lu and on the control parameter luqt^e at 
the critical value, indicates a switch between two dynam- 
ical regimes which we call a Quantum Dynamical Phase 
Transition. 

In the swapping phase the observed frequency lu is fi- 
nite. According to Eq. l(Ti)l. if cj tse ^> 1 it coincides 
with cjoj indicating a weakly perturbed evolution. As one 
approaches the critical value luqtse^^, w decreases van- 
ishing at the critical point. Beyond that value lies the 
Zeno phase where the swapping freezes (lu = 0). 

The "decoherence" rate observed from the attenuation 
of the oscillation is: 



l/r = 1/ (2r SE ) for lu q > 1/ (2r SE ) 



(16) 



The dependence of the first term in Eq. (|11[> on 1/tse 
describes the decay of the initial state into system- 
environment superpositions. This decay is instanta- 
neously compensated by the "reinjection" term which 
being "in-phase" with the Rabi oscillation ensures 
< l/ r SE- Beyond the critical value, ojqtse < 1/2, the 
decay rate in Eq. I|13|) bifurcates in two damping modes 
1/ (2tse) ± V- The slowest, for uj < 2t~b> 1S: 



i 

2t se 



1 - V 1 - ( 2w 0"T SE ) 



J T SE- (I 7 ) 



This manifests the Quantum Zeno Effect: the stronger 
the interaction with the environment, the longer the sur- 
vival of the initial state. The critical behavior of the 
observables lu and 1/r^ is shown in Fig. [21 a) and b). 

In Fig. |3 different colors label excess or defect in the 
occupation of state A with respect to equilibrium. The 
hyperbolic stripes show that the swapping period T — 
2it/lu diverges beyond a finite critical value wo^se = 1/2, 
evidencing the dynamical phase transition. Near the crit- 
ical point T ~ (T c /V5) (1 - T /T§y 1/2 , where T = j£ 
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FIG. 2: (Color online) a) Swapping frequency as a function of 
the control parameter ujotse- Departure from the dashed line 
is a consequence of the environment, b) Decoherence rate as 
a function of ljotse- At uiqtse = 1/2 there is a bifurcation 
into two decay modes. 



and Tq = 4-7TTSE is the critical natural period. The criti- 



cal exponent is —1/2. Typical dynamics illustrating both 
phases are shown at the bottom. Both start with a 
quadratic decay which is beyond the FGR and results 
from the phase 4> m Eq. (|T3|) . 

Related dynamical regimes are predicted for different 
dissipative two level systems as the spin-boson model ini- 
tially addressed by Chakravarty and Leggetl*^. While in 
that model the memory effects are fundamental, our sys- 
tem allows one to reach the strong coupling limit in a 
case where they are not relevant. Hence, both, to and r^,, 
get simple but non-analytic expressions in terms of the 
control parameter manifesting the critical aspects of the 
transition to the Zeno phase. This transition has been 
observed in NMR as an abrupt drop in the frequency 
and relaxation rate 20 ' 21 . The present theory also maps 
to these spin systems2I. 

In summary, a microscopic model for an electron in a 
two state system coupled to an environment through a 
many-body interaction shows that by sweeping wo^se be- 
low the critical value ^ the oscillatory dynamics freezes. 
While a crossover to a Zeno regime is expected for strong 
SE interaction, a non-analytic point separating a finite 
parametric region of dynamical freeze from a swapping 
phase defines a quantum dynamical phase transition. 
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FIG. 3: (Color online) Contour plot of the particle density 
function G AA (t) as a function of time and u>o in units of tse- 
The vertical line indicates the critical value where the oscilla- 
tion period diverges. Panels at the bottom show the behavior 
of G AA (t) f° r values of wotse within the Zeno phase (left 
panel) and the Swapping phase (right panel). 



